Magnetic Field Effect for Two Electrons in a Two Dimensional Random Potential 
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We study the problem of two particles with Coulomb repulsion in a two-dimensional disordered 
potential in the presence of a magnetic field. For the regime, when without interaction all states are 
well localized, it is shown that above a critical excitation energy electron pairs become delocalized 
by interaction. The transition between the localized and delocalized regimes goes in the same way 
as the metal-insulator transition at the mobility edge in the three dimensional Anderson model with 
broken time reversal symmetry. 
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I. INTRODUCTION 



The interplay of disorder and interactions in electronic 
systems is a central problem in condensed matter physics 
W . Two-dimensional (2D) systems are of a particular in- 
terest, since the scaling theory of localization || predicts 
that noninteracting electrons are always localized in a 
2D disordered potential, while a metal-insulator transi- 
tion has been reported in transport measurements with 
2D electron and hole gases ||. The study of this many- 
body problem is very complicated, both for analytical 
and numerical analysis. It is therefore highly desirable to 
have some relatively simple models which could be solved 
and would lead to a better understanding of the effects 
of interactions in the presence of disorder. The problem 
of two interacting particles (TIP) in a random potential 
has received much attention in the last few years. It has 
been shown [Q that the TIP can propagate coherently 
on a length l c which is much larger than the one-particle 
localization length l\ that can lead to an enhancement of 
transport || . This problem has been studied recently by 
different groups in one P|-|lCfl and two dimension 1 1 - 13|] 
and it has been understood that the pair derealization is 
related to the enhancement of interaction in systems with 
complex, chaotic eigenstates. The delocalization factor is 
determined by the density of two-particle states coupled 
by the interaction, pi, and by the interaction induced 
transition rate, r e . At n e — T e p2 ~ 1 the interaction ma- 
trix elements become comparable with the two-particle 
level spacing and the collisions between particles give a 
strong increase of the ratio l c jl\. 

Most studies of the TIP problem have considered a 
short range interaction. In this case, when two par- 
ticles are localized at a distance R 3> h, the overlap 
of their wavefunctions is exponentially small, and such 
states are localized much in the same way as in the non- 
interacting case. On the other hand, when the average 
distance between particles is larger than h, the screen- 



ing of charges in the problem with a given charge den- 
sity is problematic. Due to that for such a regime it is 
natural to consider the bare Coulomb interaction in the 
simple TIP problem. Recently, it has been shown E3] 
that the Coulomb repulsion can delocalize two particles 
(electrons) in a two-dimensional disordered lattice, even 
if the particles are separated by a distance R » l\. The 
delocalization of two-electron states takes place in a way 
similar to the single particle Anderson transition in three 
dimensions. Indeed, the pair center of mass moves in the 
2D plane, while electrons rotate around it that gives an 
effective third dimension. The rotation goes on a ring of 

4/3 

width li and radius R oc // fixed by energy conserva- 
tion. As a result the two-particle states are delocalized 
for excitation energies e > e c oc l ± 4 ^ 3 (n e > 1) [fl3|| . This 
expectations have been confirmed numerically |13|| by a 
study of the level spacing statistics P(s), which displays 
a transition from the Poisson distribution (for e < e c ) to 
the Wigner-Dyson distribution (for e > e c ). It was also 
found that at a critical point (e = e c ) P(s) statistics is 
close to the distribution found in the 3D Anderson model 
at the mobility edge fl4|,|l5| . 

In this paper we consider, for the first time, the effect 
of a magnetic field on the TIP problem with Coulomb 
repulsion in two dimensions. We summarize our findings 
as follows: (1) wc numerically compute TIP wavefunc- 
tions and give a direct evidence that the Coulomb inter- 
action leads to the delocalization of excited states; (2) we 
show that with the increase of the excitation energy the 
level spacing statistics P(s) exhibits a transition from 
the Poisson distribution to the Wigner-Dyson distribu- 
tion and that at the critical point P(s) is similar to the 
critical statistics found in the 3D Anderson model with 
broken time reversal symmetry |lq ]. 

The paper is organized as follows. The model is intro- 
duced in Section [nj. In Section III we review the analyti- 
cal arguments developed |l3| for the TIP problem in two 
dimensions with Coulomb repulsion and we discuss the 
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influence of a magnetic field on this theory. In Section 
iy|and[y| we discuss our numerical data for this problem 
when the time reversal symmetry is broken by a mag- 
netic field. A number of typical examples of interaction 



induced pair derealization is shown in Section IV. The 



transition in the level spacing statistics from the Poisson 
distribution to the Wigner-Dyson distribution is analysed 
in Section [v|. There we present the comparison of our re- 
sults with the data for the 3D Anderson transition with 
time reversal symmetry broken by a magnetic field. In 



Section VI we present the summary of the results. 



II. THE MODEL 

We consider two particles with Coulomb repulsion in 
a two-dimensional disordered square lattice, in the pres- 
ence of a constant magnetic field perpendicular to the 
plane. We restrict our investigations to the triplet case, 
which corresponds to the study of two spinless fermions. 
The singlet case, investigated in one and two dimensions 
for the on-site Hubbard interaction , should give sim- 
ilar results. The Hamiltonian of the model reads: 



Hint 



(1) 



<r,r'> 



The vectors r = (x, y) denote the L x L sites of a 
square lattice with periodic boundary conditions applied 
in both directions, c\ (c r ) creates (destroys) an elec- 
tron in the site r. The occupation number at the site 
r is n r = c\c r . The uncorrelated random energies E r 
are distributed with constant probability within the in- 
terval [—W/2,W/2], where W denotes the magnitude 
of the disorder. The nearest neighbors hopping terms 
on the square lattice include the magnetic field being 
V TT i = V"exp(±i27ray) for r — r' = (±1,0), while for 
r r' = (0, ±1) they are V T: r' = V. This choice cor- 
responds to the Landau gauge for the vector potential, 
A = (— Bya, 0, 0), with the magnetic field B perpendic- 
ular to the plane. The number of flux quanta per unit 
cell of the lattice is a = eBa 2 /h — Ba 2 /(f>o and in the 
following the lattice spacing constant a is taken to be 
unity. The magnetic field is chosen to be commensurate 
with the lattice, i.e. a = k/L, with k integer. The last 
term in (Q) gives the interaction: 



Hi, 



U 



V 



n r n r 



2 t", |r - 



(2) 



where U is the strength of the Coulomb repulsion and 
|r — r'| is the interparticle nearest distance computed on 
a 2D torus. 



III. ANALYTICAL ESTIMATES 

We consider the case with the average distance between 
electrons R = \r\ — r 2 | much larger than their one-particle 
localization length: R ^> l\. In the localized regime the 
one-body Anderson localized orbitals can be represented 
in the lattice basis as: 



(p a (r) w — exp 



i8 a (r) 



(3) 



where r a marks the center of the localized a-th single 
particle eigenstate, and 9 a (r) is a random phase. Interac- 
tion matrix elements between noninteracting two-particle 
eigenstates \a(3) and |7<5) are given by 



( 1 S\H mt \aP)=Ql 5 -Q% 



(4) 



with 



nl5 _ T .^ ^(r)^(rQ^(r)^(rQ 

Due to one-particle exponential localization, Coulomb 
repulsion can induce electron jumps only inside the lo- 
calization domain of size l\. Therefore, when R 3> h, 
it is possible to expand the interaction for electron dis- 
placements Ari, AV2 of typical length l\ near their initial 
positions i"i,r2. The terms up to the first order in the 
expansion of the Coulomb potential give only mean-field 
corrections to the one-particle potential. The first term 
beyond mean field has a dipole-dipole form, and is of the 
order of: 



^~-^A ri -Ar 2 ~^. 



(6) 



This gives dipole-dipole matrix elements between nonin- 
teracting eigenstates: 

(QlfiU ~ - J E Ari • ^M^M^inWs^). 

ri ,r2 

(7) 

The sum in (0) runs over l\ sites for each electron, so that 
in total the sum contains of the order of if terms with 
random signs. Each term is of the order of Ifcj) 4 ~ l^ 2 . 
As a result, the typical dipole-dipole transition matrix el- 
ement in the ergodic approximation and with eigenstates 
given by (||) is of the order of 

(8) 



->typ 

!dd 



On the basis of this result we can estimate the typi- 
cal interaction induced transition rate T e between nonin- 
teracting two-particle eigenstates by means of the Fermi 
golden rule: 



U 2 lj 
R 6 V 



(Q 



dd 



P2 



(9) 
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Here we took the density of states coupled by interaction, 
in the middle of the energy band of width B ~ V, be- 
ing p 2 ~ lf/V. Indeed, due to localization, one-electron 
jumps on a distance larger than fa give exponentially 
small matrix elements and these transitions can be ex- 
cluded from consideration. The mixing of two-electron 
states takes place when 

/ 7T/4 \ 2 

«' = T °P»~(vW) >!> ( 10 ) 

that corresponds to R < fa^fa/V) 1 ^ 3 . For U ~ V one 
gets R < l^ 3 , and the condition R 3> fa is still satisfied 
when fa 3> 1 [0. Therefore the physical picture is qual- 
itatively different from the case of short-range screened 
interaction, where mixing is possible only for states at a 
distance R < fa. For n e > 1 the pair jumps on a typical 
length fa and its diffusion rate is 

D e ~ Z 2 r e ~ -J^. (11) 
n 

The transition from localization to pair diffusion takes 
place in a way qualitatively similar to the Anderson 
model in 3D. Indeed, the pair center of mass can move in 
the 2D plane and in addition the electrons diffusively ro- 
tate around it in a ring of radius R and width fa, keeping 
their Coulomb energy E ee ~ U/R constant. The num- 
ber of effective sites in the third direction, M e f w nR/fa, 
is given by the number of circles of size fa in the ring. 
Therefore, following standard results for the quasi-2D 
Anderson model jL8|, the pair localization length l c is 
given by 

^ ~ exp(M e/52 ) ~ exp f —^j , (12) 

where 172 ~ K e is the two-particle conductance |^| and the 
above estimate is valid in the metallic phase for the cor- 
responding 3D Anderson model (n e > 1). Since R ~ l^ 3 
when K e ~ 1 (for U ~ V"), at the transition the TIP lo- 
calization length jumps from l c ~ fa to an exponentially 
large value 

l c ~ fa exp(wl\^ 3 ). (13) 

The TIP diffusion will be eventually localized due to the 
finite number of planes in the third direction. However, if 
disorder is not too strong (fa 3> 1), Coulomb interaction 
gives raise to an exponentially sharp localization length 
enhancement, with a "critical" behavior similar to the 3D 
Anderson model for finite system sizes l c >L»ii. We 
remark that, due to disorder, the ring is not rigid and the 
electronic motion adapts to the fluctuations of the ran- 
dom potential so that the total energy remains constant 
even with some variation of the ring radius. This effect 
should increase the number of planes M e / of the effective 



quasi-2D Anderson model, giving a stronger derealiza- 
tion effect. 

Since the excitation energy e is related to the pair 
distance, e ~ U / R, the condition of pair derealization 
(k e > 1) implies 

e > e c oc i; 4/3 ; (14) 

the scaling relation e c l\^ = const is in agreement with 
the numerical results obtained in Ref. |[13i . 
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FIG. 1. Single particle density of states pi (top) and the 
inverse participation ratio £1 ~ if (bottom) in the 2D Ander- 
son model as a function of energy, for system size L — 24, 
disorder strength W = TV, rescaled magnetic field a = 1/6 
(full line) and a = (dashed line). Data are averaged over 
Nr = 10 3 disorder realizations. 

Finally we discuss the influence of a magnetic field on 
the theory presented in this Section. For a typical mag- 
netic field corresponding to a — 1/6 flux quanta per pla- 
quette, the disorder strengths W — TV or W — XGV is 
strong enough to mix different Landau levels, as it is il- 
lustrated in Fig.|l| top: the single-particle density of states 
at W = TV is only slightly changed with respect to the 
zero magnetic field case. In Fig.[j] bottom we evaluate the 
energy dependence of the single particle inverse partici- 
pation ratio £1 (IPR): 

6 = i/£k>(r)r\ h-VTi- (15) 

r 

This convenient characteristic £1 determines how many 
sites contribute to an eigenstate and is simply related 
to the localization length fa. The magnetic field gives 
an increase of £1 in the middle of the energy band. In- 
deed, it is known from weak localization p8[ that, in 
the presence of a magnetic field, coherent time-reversed 
paths are eliminated and therefore backscattering is sup- 
pressed. On the contrary, the magnetic field shrinks the 
band and at W — TV there is a slight reduction of the 
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density of states near the band edges, which brings about 
a decrease in l\ |l9|. For the parameters of Fig.[l] the 
results found in fll3| ] give the critical derealization en- 
ergy e c w 1.2V, counted from the energy of the ground 
state (Fig.lc there), that corresponds to E w — 3.8V in 
Fig.l. At this energy = 1/6) sw = 0) (see 

Fig.l bottom) that implies ii(a = 1/6) « Zi(ev = 0) and 
hence we expect that the critical energy for delocalisa- 
tion transition in a presence of a magnetic field remains 
approximately the same: e c (a = 1/6) «s e c (a = 0) pOfl. 



lattice basis. ^From this we extract the one-particle prob- 
ability, 



(16) 



and the probability of interparticle distance, 

/ dfe (R)=^F fe (R + r 2 ,r 2 ), (17) 




FIG. 2. Probability distributions / (left) and f d (right) 
from Eqs. (16) and (17) for two interacting particles in a 
2D lattice of size L = 30, with disorder strength W = 10V 
and rescaled magnetic field a = 1/6. Top: Coulomb re- 
pulsion U — 2V , ground state, number of sites occupied 
by one particle given by IPR is £ s = 8; middle: U = 2V , 
rescaled one-electron excitation energy e/B = 0.71 (B — 4V), 
£ s = 320; bottom: same excitation energy but at U — 0, 

= 14.7. 



with R = ri — r 2 . 

Typical examples of probability distributions are 
shown if Fig.||, at W — 10V for a system size L = 30. 
They clearly show that the two-particle ground state 
(Fig.^j top) remains localized in the presence of interac- 
tion, with the particles sitting far from each other in or- 
der to minimize Coulomb repulsion. Similar conclusions 
apply to low-energy eigenstates. On the contrary, for 
higher excitation energies (e/B = 0.71 in Fig.|| middle, 
where e = SE/2, the excitation energy of the TIP eigen- 
state 8E is counted from the ground state and B = AV 
is the band width at W = 0) the probability distribution 
/ spreads over the whole lattice, while fd shows a hole 
at small R and a depletion for large R. The first prop- 
erty is a simple consequence of Coulomb repulsion, while 
the second one is in agreement with the general discus- 
sion of the model (0). Following the analytical arguments 



|13| summarized in Section III, we believe that this ring 
structure (see also Fig.||top) would become more evident 
at larger system sizes, with maximum interparticle dis- 
tance R m ax ~ L/2 ^> l^ 3 ||l[] . It is impossible to fully 
satisfy such a condition within the system sizes numeri- 
cally tractable, if one considers that the condition l± 3> 1 
should be satisfied at the same time [p2| . We stress that 
this pair derealization takes place in a regime of strong 
localization for the one-particle wavefunctions. This is 
demonstrated in Fig.^J bottom and Fig.|| bottom, which 
show the probability distributions for the noninteracting 
problem (U = 0) at the same excitation energy. As a 
quantitative measure of the interaction induced charge 
derealization one can take the inverse participation ra- 
tio £ s for the one-particle probability /: 



& = l/£/ 2 (r0. 



(18) 



IV. TIP DELOCALIZATION 



In order to study the eigenstate properties of our model 
with interaction, we diagonalize numerically the Hamil- 
tonian (0). In this way we determine the two-particle 
probability distribution ^(ri,^) = |^fc(ri, r 2 )| 2 for the 
k-th eigenfunction \I>fc(ri,r 2 ) = (rir 2 |^ fc ) written in the 



In this way £ s gives the number of lattice sites occu- 
pied by one particle in an eigenstate. For the case of 
Fig.^, Coulomb interaction does not significantly change 
£ s Ri 8 for the ground state, while for an excitation en- 
ergy e/B = 0.71 there is a huge derealization effect from 
£ s = 14.7 at U = to £ s = 320 at U = 2V. 
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FIG. 3. Density plots for the case of Fig.^| middle corre- 
sponding to two top plots (U = 2V) and for the case of Figj^ 
bottom corresponding to two bottom plots (U = 0). Blue 
corresponds to the minimum of the probability distribution 
ad red to the maximum. 
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FIG. 4. Dependence of rju on the rescaled one-electron en- 
ergy e/B for W = 7V, a = 1/6, U = 2V, system sizes L = 6 
(circles, number of disorder realizations Nr = 2 x 10 4 , number 
of spacings for each point of the graph Ns > 4 x 10 ), L = 12 
(stars, N R = 5 x 10 2 , iV s > 1.1 x 10 4 ), L = 18 (diamonds, 
Nr = 10 2 , N s > 1.4 x 10 4 ), and L = 24 (triangles, N R = 10, 
N s > 3 x 10 3 ). 



V. SPECTRAL STATISTICS 

The qualitative change of the structure of the eigen- 
states also leads to a change in the level spacing statistics. 



In the one particle problem spectral fluctuations proved 
to be a very useful tool to characterize the 3D Anderson 
transition |23|. Localized wavefunctions yield uncorre- 
lated spectra with Poisson statistics, characterized by a 
distribution P(s) of the energy spacings between succes- 
sive levels going to 



P P (s) = exp(-s) 



(19) 



when the system size L ^> l\. Delocalized wavefunctions 
yield correlated spectra and Wigner-Dyson statistics with 
the Wigner surmise 



Puis) 



32s 2 



■ exp 



As' 

7T 



(20) 



in the absence of time-reversal symmetry. The striking 
advantage of such an approach is that it deals only with 
the spectrum and does not involve heavy numerical cal- 
culations of conductivity or eigenfunctions. 




FIG. 5. Transition from the Poisson distribution (dashed 
curve) to the Wigner surmise (full curve) in the level statis- 
tics P(s) for the case of Fig]| at L = 18: < e/B < 0.1 
(diamonds, rju = 0.89, N s = 4.5 x 10 3 ) and 0.47 < e/B < 0.5 
(circles, nu = 0.06, iV s = 3.5 x 10 4 ). 

To analyze the evolution of the P(s) distributions with 
respect to the excitation energy, it is convenient to use 
the parameter 



Vu 



JZ v [P(s)-Pu( 3 )]ds 
C[P P (s)-P u (s)]ds' 



(21) 



where su — 0.5076... is the first intersection point of 
Pp(s) and Pu(s). In this way Pp(s) corresponds to 
rju = 1 and Puis) to rju = 0. The dependence of r\u on 
the one-electron excitation energy e = 5E/2 is shown in 
Figfl This figure shows that at fixed interaction U = 2V, 
disorder W — 7V, and rescaled magnetic field a = 1/6, 
curves at different system sizes 6 < L < 24 intersect at 
e c /B w 0.33, with r) Uc ~ 0.26. While for e < e c the level 
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spacing statistics approaches the Poisson limit (rju — > 1) 
when the system size increases, for e > e c the tendency 
is towards the Wigner-Dyson distribution [rju — » 0). Wc 
note that the critical excitation energies at VP = 7V and 
VF = 10V (e c « 0.67, data not shown) are similar to the 
values found in Il2], in agreement with the the expecta- 
tions of Section [II. 




FIG. 6. Level statistics P(s) at the critical point, for 
U = 2V and a = 1/6: W = TV, 0.32 < e/B < 0.35, L = 6 
(circles, 2V S = 8 x 10 4 ) and L = 18 (diamonds, TVs = 2 x 10 4 ); 
W = 10 V, 0.65 < e/B < 0.70, L = 12 (squares, JV fl = 2x 10 4 , 
TVs = 10 5 ). The thick line shows the critical P(s) in the 3D 
Anderson model in the presence of a magnetic field (data 
taken from Fig.l of Ref. [16]), with a large-s exponential 
decay P(s) oc exp( — ks), k — 1.87, shown in the inset. 

The level spacing statistics P(s) is shown in Fig.^Jnear 
the ground state and for excitation energies e > e c , for a 
system size L = 18. The transition from the Poisson to 
the Wigner-Dyson statistics is evident. 

The P(s) statistics near the critical point e = e c is 
shown in Fig|| The curves at different system sizes dis- 
play a size-independent intermediate distribution, which 
exhibits level repulsion P(s) oc s 2 at small s and a 
Poisson-like tail P(s) oc exp(— as), with a w 1.9. The 
close agreement between these distributions and the crit- 
ical statistics found in the 3D Anderson model with bro- 
ken time reversal symmetry at the mobility edge, taken 
from Ref. |lq ] , supports the analytical arguments given 
in Section |III| . This 3D critical statistics also gives a good 
approximation for the 2D critical statistics at W = 10 V 
(see FigJ^), with rju c w 0.29. The small difference in 
the t\tj c values at W — IV, IQV (already observed at 
a = 0, see Ref. Jl3|) could be attributed to the fact that 
the jump in the localization length at the "transition" , 

1 /3 

lc/h ~ exp(7rZ 1 ), is not sufficiently sharp. This is due 
to the not very large values of l± accessible for numeri- 
cal simulations. Indeed, the investigation of cases with 
larger l\ would require a significant increase of the sys- 



tem size, in order to satisfy the condition L > R » l^ 3 . 
We also note that the finite statistics and the limited sys- 
tem sizes prevent us from precisely evaluating the critical 
excitation energy e c and the critical value r\u c . 



VI. CONCLUSIONS 

We have shown that Coulomb interaction can delo- 
calize electron pairs in a 2D disordered potential in the 
presence of a magnetic field, above a critical excitation 
energy, in a way similar to the Anderson transition in 
3D. The close relation between these two transitions is re- 
flected in the close similarity of level statistics at the crit- 
ical point. The results obtained in this papers should be 
relevant for experiments at small electron density and/or 
large disorder fluctuations, when the distance between 
electrons is larger than the size of their localization length 
in the absence of interaction. 

We thank Isa Zharekeshev for the possibility to use the 
data of Ref. @ , and the IDRIS in Orsay and the CICT 
in Toulouse for access to their supercomputers. 
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tion coming from the expansion of the Coulomb repulsion 
near R ~ U/e. 

We note that L = 30 corresponds to a large Hilbert space 
dimension Nh ~ 4 x 10 s and that numerical diagonaliza- 
tion of such matrices is possible thanks to the Lanczos al- 
gorithm, with special care to find eigenvalues in the inte- 
rior of the spectrum, see J. Cullum and R.A. Willoughby, 
J. Comp. Phys. 44, 329 (1981). 
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